Quantum simulation of multiple-exciton generation in a nanocrystal by a single photon 
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We have shown theoretically that efficient multiple exciton generation (MEG) by a single photon 
can be observed in small nanocrystals (NCs). Our quantum simulations that include hundreds 
of thousands of exciton and multi-exciton states demonstrate that the complex time-dependent 
dynamics of these states in a closed electronic system yields a saturated MEG effect on a picosecond 
timescale. Including phonon relaxation confirms that efficient MEG requires the exciton-biexciton 
coupling time to be faster than exciton relaxation time. 

PACS numbers: 



Solar light would be an important souree of clean and 
renewable energy if the efhciency of inexpensive solar 
cells could be increased. Increased efficiency can be 
achieved through carrier multiplication: Photo-generated 
carriers, whose excess energy is greater than the energy 
gap, can create secondary electron-hole pairs via impact 
ionization of the filled band. Through this process, two 
(or more) electron-hole pairs are collected from each pho- 
ton instead of just one. This process is very inefficient in 
bulk semiconductors, where impact ionization has a very 
low probability and carrier thermalization, which always 
competes with impact ionization, is much faster. 

As first suggested by Nozik, impact ionization may ef- 
fectively compete with cooling in NCs, due to the en- 
hanced rate of inverse Auger processes for carrier multi- 
plication and the "phonon bottleneck" suppression of car- 
rier relaxation, leading to efficient MEG [l[. Soon after 
this publication, Schaller and Klimov [^ observed ultra- 
efficient MEG by a single photon in PbSe NCs, using 
band-edge transient absorption measurements. Later, ef- 
ficient MEG was observed by many groups using different 
techniques in NCs of man y se miconductors: PbSe [3H6|, 
PbS [113, Si i^CdSe (i^lll, InAs [Hill, and in car- 
bon nanotubes |15| . At the same time, some groups were 



u 



not able to observe MEG in CdSe [16l and InAs |14l.ll7| 
NCs and found that the efficiency of MEG measured in 
PbSe NCs p was appreciably smaller than that reported 
earlier. The diverse experimental data on the MEG effi- 
ciency are now converging to more modest values for PbS 
and PbSe NCs O], but MEG in NCs has been shown to 
be significantly more efficient than impact ionization in 



bulk semiconductors [19|, l20| . 



Previous attempts at a theoretical understanding of 
the enhanced MEG provide only estimations of the 
MEG efficiency observed in NCs 0, ^^^. A self- 
consistent theory of this phenomena requires a currently 
non-existent theoretical description of both the relax- 
ation mechanisms for and couplings between the highly 
excited exciton and multi-exciton states in NCs [2J] . To 
explain the high efficiency of MEG in NCs, the coherent 



superposition model , based on the strong quasi-resonant 
coupling between exciton and multi-exciton states in a 
NC, was proposed Jsl |2J|. Non-coherent models for effi- 
cient MEG in NCs |2ll423ll26j are based on the important 
observation that the density of biexiton states is signif- 
icantly larger than the density of exciton states at the 
same energy [21| , and the density of trion states is partic- 
ularly important for efficient MEG [25| . The calculations 
of MEG efficiency in all of these non-coherent models are 
based on Fermi's Golden Rule, which requires the final 
biexciton state to decay much faster than the rate of the 
exciton-biexciton transition, an assumption that has not 
been justified either experimentally or theoretically. 

In this letter we unify both approaches and consider a 
single-photon excitation coherently coupled with multi- 
exciton-states in a NC within a full quantum-state evo- 
lution approach. The time-dependent dynamics of our 
modeled systems is described using a large multiple- 
configuration basis representation of the many-electron 
Hamiltonian, including energy non-conserving exciton 
and biexciton decay channels. The calculations show that 
even in a closed, energy-conserving electronic system, the 
excitation becomes predominantly multi-excitonic on a 
picosecond time scale. The initial single-photon excita- 
tion is dispersed into the dense multi-exciton state space 
of the NC. 

We consider a single-photon excitation of a spheri- 
cal PbSe NC, where efficient MEG has been reported 
l2-|6l- The energy spectrum of electrons and holes at 
the band edges of four equivalent L valleys of bulk 
PbSe is described by a four band k ■ p model j27| . 
In PbSe NCs each electron, n^L^m , ^^i^ hole, n^L^m j 
state within the four-band effective- mass model [281, is 
characterized by the spatial angular momentum L ~ 
0{S),1{P),2{D),3{F), ... , the total angular momentum 
j = L ± 1/2 , which is the sum of the angular momen- 
tum and the spin, the projection of the total momentum 
m = ±1/2, ±3/2, ... ± J , and the spatial inversion parity 
TT . The energy levels of the same symmetry were cal- 
culated using energy band parameters from Ref. y and 



labeled by a level number, n^'^ . The four envelope 
components, Fi{r) , [i = 1,2,3,4) of the correspond- 
ing single-particle wavefunctions can be expressed in the 
form: 



(a) t = 0.64 fs 



F,(r)=Cf'™'Vp„(r)r":-(^,'^), 



(1) 



where C^ '"^'^ are constants, fp{r) are specific radial 
functions expressed in terms of the spherical Bessel func- 
tions, LJ ~ L OT L + 1 and mf = m ± 1/2 (depending 
upon the values of i and tt), and Yp^{9,(j)) are the 
spherical harmonics. 

The Coulomb coupling of an electron or hole with 
multi-electron-hole excitations at the same energy, which 
are disallowed in the bulk due to momentum conserva- 



tion, plays an important role for MEG in NCs [2J|. We 
generate a Hamiltonian matrix in the basis of many- 
electron states that are Slater determinants of single- 
particle eigenstates. A many-electron state in this ap- 
proach is a set of occupied conduction states (electrons) 
and unoccupied valence states (holes) with an arbitrary 
but consistent sign convention for the Slater-determinant 
permutations. We include many-electron basis states up 
to some chosen energy cut-off that we adjust until the 
results converge (see Fig. [1]). We modeled the Coulomb 
interaction using effective dielectric constants for the NC, 



Kg , and surrounding media. 



i-g , 



leading to the direct 



Coulomb interaction between electrons at positions r 



and r' in the NC, Vcir,r') 



'/i^. 



the surface mediated Coulomb potential 



a/{Kr'\r 
-1 = 2(k, 



aV/r'^l) 

- Kg)/ Kg[Ks 



- r'\), and 
F,(r,r') = 



where a is the NC radius and 
+ 2Kg) . We also account for the 



electron and hole interaction with their image potentials 
described by T4(r) = 14(r) = 0.5V"s(r, r) . 

The diagonal elements of this Hamiltonian matrix are 
taken to be sums of the singie-particle-state energies (in- 
cluding image potentials) and the Coulomb interactions 
among electrons and holes (including the exchange inter- 
actions from Fermi permutations) . Off-diagonal elements 
couple basis states that differ in either two or four sets 
of single-state occupations via the Coulomb interaction. 

Because the single-particle states in the spherical NCs 
are represented in terms of spherical harmonics, it is most 
efficient to express the Coulomb potentials V(7(''", '''') and 
Vs{r, r') in terms of spherical harmonics [29|: 



Vs{ry) 



., ^2/ + l.^+i^ 
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4^e2^ 1 (rr''^ 
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1 Tz(0,0,0',0'), 



1=0 



where T,(0, 0, ^',0') = E™ >1™(^', '^')i"/-(^, 0) , r and 
r' are represented by {r, 0, 0} and {r',^',(/)'} in spher- 
ical coordinates, r> = max (r, r') , and r< = min (?', r') . 
Using these expressions, we have performed integrations 
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(c) t = 655.36 fs 
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FIG. 1: Time dependent evolution of the 2''P-^ 
citation created in the 2 nm radius PbSe NC witli Ks ~ Hg = 5 
by a single Ej = 2.95 eV photon calculated for the closed- 
system. The upper panels (a,b,c) show the occupation proba- 
blities of exciton, bi-exciton and tri-exciton states over 0.2 eV 
wide energy bins at three times. Panel (d) demonstrates con- 
vergence of the time dependence of the relative number of 
excitons {Nx{t)) — 1 , defined in Eq. ([Sj, achieved by adjust- 
ing the energy cut-off'. The vertical fines a, b and c show the 
times corresponding to the three upper panels. 



involving the angular components analytically by exploit- 
ing the orthonormality of the spherical harmonics. 

The remaining radial parts of the integrals 
of the Coulomb matrix elements take the form 
Jdrjdr'frir)fn^)mr')fr{r')[ir^y/ir>y+'], 

where 



\i+i 



Jdrjdr'fr{r)fr{r)mr')fr{r'){l/rr') 
fl"'{r) are the radial functions defined in Eq. ([T]) for the 
envelope eigcnfunctions that are expressed via spherical 
Bessel functions. In our computational implementation, 
we approximate these Bessel function integrals as needed 
by Monte-Carlo sampling and store the result for future 
reference. A sufficient number of samples was taken 
to reach a specified level of precision for the matrix 
elements [30J. At a 4.8 eV cut-off used to simulate 
the 2 nm radius NCs, the basis contains hundreds of 
thousands of many-electron states, and the Hainiltonian 
matrix contains tens of millions of non-zero elements. 

To incorporate the interaction involving a single pho- 
ton, we adopt the standard effective-mass approxima- 
tions that the wavelength of the photon is much larger 
than the NC size and that the contributions of the Bloch 
functions are dominant in the integral over the electron 
momentum operator [28| . In the PbSe NCs with symmet- 
ric conduction and valence bands, we assume that the 
photons are coupled exclusively to symmetric electron- 
hole pairs, in which both the electron and the hole have 
the same quantum numbers. As a result the NC-photon 
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FIG. 2: Time dependent evolution of tiie relative num- 
ber of excitons {Nx{t)) , defined in Eq. ([3|, created by a 
single-photon excitation of various optical transitions in PbSe 
NCs with radii 2 nm (a) and 3 nm (b) calculated for the 
closed-system with Hs = Hg ~ 5 . The optical transi- 



tions 2''Plj^ 



' 2''P^{^ , VD^i^^ - \^dV?^ , rF':{l - \^f: 
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1/2 ' '- -^1/2 '- -^1/2 ' -^ -^ 1/2 -^ -' 1/2 

^T^il - 2^pV^I , rcVl - i^gVI , ^''dVI ~ 2^d^JI 
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l'J/f/^2 - I'^-ffj'/a ' are abbreviated ID , IF, 2P , IG , 
2D , and IH respectively with energies shown in units of 
the effective energy gap Ei{a) . Panel (c) shows the evolu- 
tion of {Nx{t)) created by a single-photon excitation of the 
2'^Djy2 ~ 2 -Dj/2 transition in a 2nm radius PbSe NC calcu- 
lated for various combinations of the dielectric constants of 
the nanocrystal Ks and the surrounding media Kg . 



interaction is characterized by a single parameter cor- 
responding to the coupling energy. Furthermore, our 
results are obtained in the weak-coupling single-photon 
limit, where the actual coupling strength is irrelevant 
provided that it is sufficiently small. 

Let us begin from a study of the closed system evolu- 
tion. Starting from an initial single-photon-state, elec- 
tronic states are excited and evolve through Coulombic 
interactions. In each simulated experiment, the photon 
energy was chosen to match one of the optically active 
exciton states, because the NC absorption spectra are 
not affected by the coupling between excitons and multi- 
exciton states 24|. We evolve the initial state using the 
Schrodinger equation, d\'^{t))/dt = —iiJ|^(t))/fi, , where 



H is our Hamiltonian that includes all interactions con- 
sidered above. Using a difFerental equation solver [31| for 
integration, we adjusted the precision until we obtained 
convergent results. 

Figure [T] shows the time-dependent evolution of the 



2^R 



1/2 



1/2 



1 /? 

2'^P^L excitation created by a single pho- 
ton calculated for our closed-system description. The 
upper panels (a,b,c) show the occupation probablitics 
||(fc|^(t))||^ of exciton, bi-exciton and tri-exciton states 



created from the photon at three different snapshots in 
time. The initially created excitons are replaced by the 
biexcitons in approximately 0.3 ps. Formation of tri- 
excitons at the excitation energy considered, 2.95eV^ ~ 
3.73-Ei(a) where Ei{a) is the effective energy gap, is 
limited, however, in our single valley model, which al- 
lows only two electron-hole pairs to occupy the ground 
state. Figure [TJd) shows the time dependence of the ex- 
pectation value for the relative number of excitons: 



ri=l fe6K„ 



;>o 



where /C„ is the set of n -exciton basis states. In the 
small NC with small dielectric constant considered in 
Fig. m the strong Coulomb coupling rapidly creates a 
bi-exciton with 90% probability. 

We also investigate how {Nx{t)) depends on the ex- 
citation energy and the Coulomb coupling, which in 
our model is controlled by the NC size and dielec- 
tric constants. Figure HJa) compares the formation of 
multi-excitons from a single photon excitation of dif- 
ferent optical transitions in the 2 nm radius NC. The 
2'^P^ ,'2 ^ '^^^xii transition appears to provide a thresh- 
old for efficient MEG generation in PbSe NCs: while 
lower transitions saturate at around 10 fs to a low prob- 
ability of biexciton generation, the multi-exciton number 
continues to increase when the threshold is reached. This 
is because the T^P and 2^P states are the first that can 



be strongly coupled with trions of the same energy [24 1 . 
Increasing the NC size [Fig. [DJb)] and the dielectric con- 
stant [Fig. HJc)] significantly increases the time of MEG 
generation. Both effects are the direct result of the re- 
duced Coulomb coupling. 

Note that the evolution of (iV^(t)) does not show 
strong oscillations when MEG is efficient, even for a 
closed system. Oscillations are seen only in the strik- 
ing cross-over behavior in Fig. [2] (a) and (b), which is the 
result of the different coupling strengths and detunings 
between the multi-exciton states and the exciton created 
by the resonant light. The coherent oscillations connect- 
ing the exciton and biexciton states are suppressed due 
to the large number of interacting biexciton states. Once 
the excitation evolves into the dense multi-exciton state 
space, it does not appreciably revive. 

Phonon-assisted assisted relaxation of carriers, which 
was not considered above, can strongly affect the MEG 
efficiency [2J]. The evolution of {Nx{t)) , including 
phonon decay, is described by 






-1hi*)-EE^(^-i*)I^^) 



(4) 



n k£K„ 



where r„ is the relaxation time of the n -exciton state 
that for simplicity is assumed to depend only on n . As 
states decay, they are removed from the Hilbert space 
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FIG. 3: Evolution of tlie average number of decayed-to-tlie- 
ground-state excitons {Dx{t)) created by a single-plioton ex- 
citation of the l^H 



l'^H-,1^ transition of a 2 nm radius 



'1/2 ^ -"^1/2 

PbSe NC calculated for various combinations of thermaliza- 
tion times. The inset shows the distribution of inverse cou- 
pling times connecting optically active exciton states and bi- 
exciton states, which are equal to the interaction Hamiltonian 
matrix elements, divided by h . The 1/ps scale in the inset 
allows comparison of typical coupling times to thermalization 
times. The gray region of the coupling strengths is unresolved 
at our precision and excluded from the distribution. 



in our model. The intcr-valley coupling significantly in- 
creases the density of multi-exciton states; the k -exciton 
state density is increased by at least a factor of 4^'^ for 
PbSe with its 4 valleys. For this reason, MEG should be 
even more efficient in a model that includes inter-valley 
coupling. 

In summary, our calculations unambiguously demon- 
strate that highly efficient MEG can be observed in small 
NCs. The effect is enhanced by a high density of biex- 
citon states that are strongly coupled with optically cre- 
ated excitons. Fast multi-exciton thermalization acceler- 
ates the formation of multi-excitons in the ground state 
which should improve efficiency of extraction of electron 
hole pairs from NCs. 
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and populate D„ (f ) , denoting each "ground" n -exciton 
state and governed by £)„ = X^fcGK; 11(^1*) IP/'''" ^'^^"^ 
D„(0) = 0. 

To find the relative number of excitons follow- 
ing decay, we take the weighted average {Dx) = 
J2n "n-Dn/ Y.m ^m ■ We plot tliis average in Fig. [3] for 
various r„ lifetimes. The introduction of the long and 
equal relaxation times ti = r2 = ra = 1 ps has little im- 
pact upon the time dependent evolution: the system re- 
mains almost closed. Indeed, these times are much longer 
than an average time of the exciton-biexciton coupling as 
one can see in the inset of Fig. [3] Reducing the biexciton 
and trion decay times to T2 = T3 = 100 fs increases the 
rate of MEG in agreement with Ref. [2J] . It is interest- 
ing to note that the saturation value for MEG depends 
only upon ti . If this time is shorter than the exciton- 
biexciton coupling time, ~ 300 fs, MEG is significantly 
suppressed. 

Our calculations of the dynamics and efficiency of 
MEG, unfortunately, cannot be compared with experi- 
mental data measured in PbSe NCs because they were 
based on a single-valley model. The lack of a precise 
knowledge of inter-valley coupling does not allow us to 
take this effect into account. The single- valley approxi- 
mation should significantly underestimate the density of 
multiple-exciton states. For example, in real PbSe NCs 
one can have eight electron-hole pairs at the ground exci- 
ton state instead of only the two ground excitons allowed 
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